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ABSTRACT 



This paper presents the results of a comprehensive wide field search for tran- 
siting "Hot Jupiter" planets (gas giant planets with an orbital period ld<P<16d) 
in the globular cluster 47 Tucanae. Motivated by the detection of the transit in 
HD209458 and the apparent lack of planetary detections in the core of 47 Tuc by 
Gilliland and coworkers, this work further addresses the question of giant planet 
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frequency in 47 Tuc by observing, from the ground, a 52'x52' field centered on 
the cluster. Hence this work is most sensitive to the uncrowded outer regions, 
where the stellar densities are significantly lower than the core, and concentrates 
on 21,920 main sequence stars within 2.5 magnitudes of the cluster turnoff (and 
hence approaching solar in mass) Our work comprises the largest ground-based 
transit search of a globular cluster to date, incorporating a 33-night time-series 
which allows us excellent sensitivity to detect Hot Jupiter planets. 

Detailed Monte Carlo simulations incorporating the actual temporal sampling 
and photometric precision of the data predict that seven planets with an orbital 
period range of 1— 16d should be present in our dataset, if 47 Tuc has the same 
planetary frequency observed in the solar neighbourhood. A detailed search 
utilizing a matched filter algorithm, developed specifically for this project found 
no transit events. This 3.3o" result is consistent with the HST cluster core null 
detection of Gilliland and coworkers. Our result indicates that system metallicity 
rather than crowding is the dominant effect inhibiting Hot Jupiter formation 
in this environment. The 33-night dataset used for this result also led to the 
detection of 100 variable stars, including 69 new discoveries, which have been 
presented in a companion paper. 

Subject headings: globular clusters: individual (NGC 104, 47 Tucanae) - plane- 
tary systems - techniques: photometric 

1. Introduction 

Using transit photometry to search for close-in extrasolar giant planets (Hot Jupiters) 
has recently yielded success. A planet of Jupiter's radius orbiting at <0.1AU will cause a 
1-2% dimming of ~2 hours duration in a solar-type main sequence parent star, if its orbital 
inclination is sufficiently edge-on as viewed from Earth. Transits allow the determination 
of a number of important system characteristics, including the orbital period, an estimate 
for the orbital inclination and the planetary radius. When coupled with radial velocity 
observations, an accurate measure of the planetary mass is determined. Hence, such transit- 
determined parameters are important for theoretical studies of extrasolar planet formation 
and evolution. 

The availability of moderate aperture telescopes over long contiguous observing windows 
has allowed many transit candidates to be discovered in massive campaigns. HD209458b was 
the first planet to be observed with the transit method (Charbonneau et al. 2000; Henry et 
al. 2000), being first discovered using radial velocity techniques (Mazeh et al. 2000). This 
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Hot Jupiter was later followed up by Brown et al. (2001) using HST/STIS, who achieved a 
photometric precision of 0.1 mmag and a temporal resolution of 80s. Using this dataset, the 
planetary radius was estimated to be Rp=1.35±0.06Rjup. With knowledge of the planetary 
mass (Mp=0.69±0.05Mj„p), an estimate was made for the density (0.35 g cm^'^). This planet 
seems to have an unusually large radius for it's mass. BarafTc ct al. (2003) discusses the need 
for an additional source of heat to account for these radius measurements, whereas Burrows 
et al. (2000) suggests that the cause is due to high residual entropy caused by the planet's 
formation very close to the host star. Recent transit discoveries have suggested that the 
extended atmosphere of HD209458b is somewhat unusual (Burrows et al. 2004), and is likely 
caused by the increased stellar flux retarding the contraction of an extended planetary radius 
from the time of formation. 

Using the 1.3m Warsaw telescope in Chile, the OGLEIII group originally discovered 46 
transit candidate systems towards the Galactic Bulge (Udalski et al. 2002). The transiting 
signals in these systems could be caused by planets, brown dwarfs, stellar companions or 
blended eclipsing binaries. A total of 137 candidates have been presented more recently by 
Udalski et al. (2003), and are the subject of vigorous follow-up programs to determine their 
likely nature. From these candidates, OGLE-TR-56b became the first planet discovered 
using the transit method (Konacki et al 2003) and confirmed with radial velocity data that 
measured the planetary mass (Torres, Konacki, Sasselov, & Jha 2004; Konacki et al. 2003). 
The estimates for the planetary parameters of OGLE-TR-56b are Mp=1.45±0.23Mjuj, and 
Rp=1.23±0.16Rjuj„ which produce a density of 1.0 g cm~^, significantly higher than that of 
HD209458b. Very recently, two new transiting planets have been confirmed by Bouchy et al. 
(2004), 0GLE-TR-113b and OGLE-TR-132b. Of these, 0GLE-TR-113b was independently 
verified by Konacki ct al. (2004). Both of these new discoveries have very short orbital periods 
(~1.2d), much shorter than the apparent 3d cutoff seen in the distribution of radial velocity 
Hot Jupiter discoveries ^, cither indicating the presence of a new class of close-in 'Very Hot 
Jupiters', or revealing a selection effect. Either way, these objects present interesting targets 
for theoretical work. Very recently, Pont et al. (2004) announced the planetary nature of 
OGLE-TR-111, and Alonso et al. (2004) announced TrES-1, the first planet transit detected 
via the TrES multi-site campaign. 

Many groups are searching for transiting planets. For example Mochejska et al (2003) 
have undertaken the Planets in Stellar Clusters Extensive Search (PISCES) targeting selected 
open clusters, and Hidas et al. (2003) have undertaken a wide field search for transits in the 
galactic plane, and expect about seven Hot Jupiter discoveries a year. Further information 
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on the status and prospects of current selected transit searches can be found in the review 
by Home (2003). 

Globular clusters are an excellent laboratory in which to search for transiting planets. 
47 Tuc is one of the brightest globulars in the sky, and is well-placed for observation from 
the Southern Hemisphere. Additionally, it is among the closest globulars (around 4.7Kpc 

(Harris 1996)), and thus well resolved by small telescopes, and contains bright main sequence 
objects for study (V=17.0+). With a suitably wide field instrument it is possible to observe 
the whole of the cluster in a single exposure, capturing tens of thousands of target stars in 
a single field of view. This maximises the temporal coverage, increasing the probability of a 
successful detection, and raises the statistical significance of any null result. The cluster has 
excellent observability at its respective season from the south, and is located far from the 
path of the Moon. 

Motivated by the detection of the transit in HD209458 and by the apparent lack of 
planetary detections in the inner, crowded region of the globular cluster 47 Tucanae by 
GiUiland et al. (2000), we further address the question of which conditions are necessary 
for giant planet frequency by observing, from the ground, a very wide field centered on 
47 Tuc. We concentrate on transit detection involving main sequence stars close to the 
cluster turnoff (MSTO), and hence on parent stars approaching solar mass. This work 
complements that of Gilliland et al. (2000), who used WFPC2 on HST to search for transits 
on an ensemble of ~34,000 main sequence stars in a field centered on the core of 47 Tuc, 
by effectively performing a "control experiment" outside the core. We search for the same 
transit signatures in the uncrowded outer regions of the cluster, utilising a field of view 
~250 times larger, to examine whether the paucity of Hot Jupiters in the core is due solely 
to the density of the stellar environment. Ours is the largest ground-based transit search of 
a globular cluster to date. We compare our results to the frequency of known Hot Jupiters 
in the solar neighbourhood. 

The frequency of planetary systems appears to be dependent on the metallicity of the 
host star (Gonzalez 1997; Laughlin 2000; Fischer & Valenti 2003), with stars of low metaUicity 
having a lower probability of planetary formation. Globular clusters allow a direct test to be 
made on this dependence. These clusters are among the oldest of astronomical objects, and 
their stars are typically therefore of very low metallicity. Our target, while metal-rich as far as 
globulars are concerned (47 Tuc [Fe/H]=-0.76 (Harris 1996)), is still metal poor compared to 
the Solar Neighbourhood. Furthermore, globular clusters allow the effect of stellar crowding 
on planetary formation to be studied. Such crowding may disrupt protoplanetary disks and 
harden typical planetary systems. Stellar encounters would be more commonplace in the 
cluster core, sampled by the GiUiland et al. (2000) HST survey, but much less common 
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in the uncrowded outer regions where this experiment is most sensitive. By considering 
the timescale needed for planet formation and migration, which is generally regarded to be 
short (<10^yrs)(Boss 2001), along with the relaxation time, which increases with increasing 
radius from the core, the region chosen for a transit search could have a large bearing on the 
probability of a detection. Short crossing times would have a larger effect on circumstellar 
material at early times and could inhibit planet formation. We test this by sampling the 
whole of the visible cluster face (except the inner 6') out to ~60% of the cluster tidal radius 
of 45'. 9 (Leon, Meylan, & Combes 2000) in one exposure. 

Section 2 describes the observations and data reduction, along with notes on the method 
used to produce the photometry and astrometry. Section 3 details the expectations for our 
transit search if the planet frequency in 47 Tuc is similar to the solar neighbourhood, with 
a brief description of the algorithm developed to carry out the search, which is presented in 
full in a companion paper. Section 4 describes the results of our search and its implications. 
Section 5 contains the summary and conclusions. 

2. Observations and Data Reduction 

Full details of the observations, data reduction, photometry and astrometry have been 
given in Weldrake et al. (2004a), hereafter WSBFa, a companion paper presenting the vari- 
able star catalogue produced from the same dataset, and so these aspects are only sum- 
marised here. Using the Australian National University (ANU) 40-inch (Im) telescope at 
Siding Spring Observatory, 47 Tucanae was observed for a total of 33 nights from 2002 Au- 
gust 22 to 2002 September 24. Approximately 80% of the time was useful for the transit 
search, with mean seeing of 2.2 arcsec. The Wide Field Imager (WFI) was used, consisting 
of a 4x2 array of 2048x4096 pixel back-illuminated CCDs, and providing a field of view 
52' on a side, with a 0.38" pixel detector scale. This field of view allows the whole of the 
visible extent of the cluster (out to 60% of the cluster tidal radius) to be imaged in a single 
exposure. The field layout can be seen in Fig.l of WSBFa. 

In order to achieve a significant increase in S/N for any given exposure time while 
keeping atmospheric dispersion to an undetectable level, a single broadband filter covering 
Cousins V and R was used. Temporal resolution was maximised as much as possible with 
300s exposures taken on average every six minutes for around 10 hours per night. The 
dataset in total comprises 1220 images of the same field centered on 47 Tuc, which have 
been used to produce time-series lightcurvcs. Image offsets were minimised at the telescope, 
so that the database of stars is maximised and the effects of CCD gaps reduced. Any images 
containing satellite trails or taken in bad conditions were discarded at the telescope. In total 
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our database contains 109,866 stellar lightcurves with a 33-night temporal extent.. 

Image reduction was undertaken using the MSCRED package within IRAF^. Standard 
reduction practices were employed, including region trimming and overscan correction, bias 
correction, flat-fielding and dark current subtraction. 

2.1. Photometry 

We produced the time-series lightcurves using a Difi^erence Image Analysis (DIA) method 
originally described as an optimal PSF matching algorithm by Alard & Lupton (1998), and 
later modified by Wozniak (2000) for use in detecting microlensing events. Rather than 
producing the total fiux per image, this method measures the difference in fiux of stellar 
objects between individual dataset frames, and a template frame. Lightcurves are generated 
in a linear fiux unit, from which the constant flux of that object on the template frame has 
been subtracted using a specified PSF. This is made to be equal across the whole dataset, 
reducing systematic effects. We estimated the template flux of each star that was visible in 
our best seeing image (FWHM=0.8"), and then converted the linear flux unit to differential 
magnitude units for the main analysis. 

The details of our photometric procedure and uncertainties are given in WSBFa. Of 
importance here is that the fleld contains a large change in crowding between the four 
inner and four outer CCDs, with a corresponding change in photometric uncertainty at a 
given magnitude limit. We have sufficient photometric quality to detect a 1.5-2.5% dip, 
typical of transiting Hot Jupiter planets, down to V=18.5 and V=19.0 in crowded inner 
and uncrowded outer regions of 47 Tuc, respectively. Our rms uncertainty is 0.04 mag for 
V=19.5 stars, allowing us to probe at least the brightest two magnitudes of the cluster main 
sequence for any orbiting giant planets. The ability of our algorithm to detect transits even 
in the presence of substantial noise has a direct bearing on the faint magnitude limit of our 
search, as discussed in the following section. 

3. Transit Expectations 

For this work, we have defined a Hot Jupiter as having an orbital period less than 16d, 
similar to that of 14d used by Basri (2004), and roughly corresponding to an apparent gap 
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in the orbital period distribution of exoplanets^. This definition is somewhat arbitrary; our 
definition comes from the hmit of orbital period that can be theoretically detected in our 
dataset. 

In order to search for specific features in the photometric time-series, maximum sensi- 
tivity can be obtained if the parameters of the features can be anticipated, and a matched 
filter algorithm used. Here wc will describe the calculation of the expected depths and du- 
rations of transits in our dataset. Ignoring the effects of stellar limb darkening, the transit 
depth can be determined from a simple ratio of the planetary and stellar radii, as derived 
by considering the maximum fractional change in the observed flux: 



where Rp is the radius of the planet, and i?* is the radius of the stellar host. Estimates can 
be made for the mass and radius of a star in 47 Tuc at any given V and V-I from Bergbusch 
& VandenBerg (2001), and are tabulated in Tablel for stars of 17.0<V<19.5 (corresponding 
My from 3.6 to 6.1), assuming a 47 Tuc distance modulus of 13.4. 

As 47 Tuc is a very old object and with an [Fe/H] of -0.7, the stellar mass at the Main 
Sequence Turn Off (MSTO) is estimated to be -O.QMq (Bergbusch & VandenBerg 2001), 
corresponding to luminosities ~2.5Lq and radius 1.5©. It therefore follows that the stars 
sampled for transits in this paper will have a maximum mass of this limit (O.QM©). Stars 
beyond the cluster turnoff (and not on the red giant branch) are physically larger and would 
therefore exhibit lower amplitude transits. Stars on the main sequence of 47 Tuc have smaller 
radii than solar-mass stars in the solar neighbourhood, leading to a deeper transit signature, 
which makes them easier to detect. 

The MSTO for 47 Tuc is at V~17.2, taken from the work of Hesser et al. (1987), and 
consistent with that found by WSBFa from their Colour-Magnitude diagram, reproduced in 
Fig.l. We take this as the bright stellar limit in our search of 47 Tuc stars for Hot Jupiter 
companions. The lower limit depends on the ability of our detection algorithm to detect 
the deepest transits that could conceivably be caused by a planetary companion. Prom 
Monte Carlo simulations involving modelled transits of varying depth, we found that our 
algorithm can detect 3% depth (0.03 mag) transits in V=19.5 stars, which typically have 
an rms photometric uncertainty of 0.04 mag, forming the lower photometric limit to our 
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search. A total of 21,920 stars lie between these bounds; these are the stars we search for 
Hot Jupiters. 

In order to determine the expected duration of a planet transit as a function of V 
magnitude along the main sequence of 47 Tuc, we must consider the relationship between 
the length of transit (rtran) and stellar parameters (see, for example Gilliland et al. (2000): 



Ttran = IAUM'^/'R^P'J,' (2) 

where and i?* are the stellar mass and radius in solar units where Porb is the orbital 
period of the planet in days, and Tfran is also given in days. This relationship includes a 
7r/4 reduction in transit duration, and hence is the expected duration for a planet which 
crosses the 'average' chord length of the stellar disk. A centrally-crossing transit has a 
duration 1.273 times longer and a grazing transit duration is 0.66 times shorter. Using Eq.l 
and 2 and the results given in Table 1, we now calculate the transits depths and durations 
(assuming the 'average' chord length) expected for any Hot Jupiters in our dataset. These 
are shown schematically in Fig.l, where we indicate the transit depth of a 1.3 Rj planet, 
with an assumed orbital period of 3.3d, (the most typical orbit for such a planet in the Solar 
Neighbourhood^) as a function of V magnitude. The upper and lower limits of the transit 
search range are defined by the location of the MSTO and the point where our photometric 
uncertainty is poorer than 0.04 mag (rms). 

Hot Jupiter transit depths are expected to vary from 0.007 mag to 0.031 mag over our 
search range, with the shallower transits having a longer duration than the deeper, due to 
the larger radii of the parent stars. These two factors have a role to play in determining the 
visibihty of a transit in the dataset. It is clear that a transit can only occur if the orbital 
plane is inclined at a suitably critical viewing angle. This inclination i must therefore satisfy: 



a cos i — + Rp (3) 

where a is the planetary semi-major axis, i is the orbital inclination, R^ is the stellar radius, 
and Rp is the radius of the orbiting planet. However, cos i can take on any random number 
from to 1, so to derive a transit probability, Ptran for a large sample of randomly oriented 
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targets, this must be taken into account by integrating over all possible orientations. We 
take Ptran to be given by Sackett (1999): 



P.„„= ^°"7';°''':-'^ =^~0.8n^) (4) 

/q ci(cosi) a « 

We have taken the upper limit to be appropriate for a grazing transit, that is, we have 
assumed that the limit of detection is defined when the outer edge of the planetary disk just 

grazes the edge of the stellar disk. The factor of 0.87 is derived using an assumed planetary 
radius of l.SRj. Using this relationship, the transit probability of a planet of 1.3Rj orbiting 
a solar radius primary, typical of the stellar radius at the centre of our search range is plotted 
as a function of orbital period (P = l — lGd) and is presented in Fig.2. The transit probability 
drops dramatically as period increases: at 2d orbital period the probability is 13.0%, whereas 
at 8d, the probability drops to 5.2%. An orbital period of 16d has been chosen as the upper 
limit of the search as any planet at this period is at the theoretical limit of displaying three 
detectable (and hence confirmed) transits within our 33-night observing window. It is of 
course possible to have only one or two transits visible in the data. However, in oTir search 
these are wholly caused by systematic effects as they all occurred at the same times in the 
dataset and were much more common on stars in the most crowded regions of the data with 
common periodicities. 



3.1. Method for transit detection. 

We employed a Matched Filter Algorithm to search for the periodic transit signature 
in our time series lightcurves. Full details of the method used are given in our second 
companion paper, Wcldrakc & Sackett (2004b), hereafter WS. The matched filter method 
was first suggested for use in transit detection by Jenkins, Doyle, & Cullers (1996), has been 
described as the best method for transit searches in the literature (Tingley 2003a,b), and has 
been used as the method of choice for several transit searches in recent years (Gilliland et al. 
2000; Bruntt et al. 2003). Our MFA was designed and produced specifically for this project 
and has been extensively tested with Monte Carlo simulations to test for and maximize the 
recoverability of modelled transits while keeping the false detection rate to an acceptable 
level. These tests involve the addition of false transits of various depths, durations and 
periods onto the actual time-series dataset. Thus, the actual time sampling and photometric 
uncertainty of the data is taken into account by our Monte Carlo recoverability simulations. 
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The method as implemented assumes a simple square-well model, which is a valid as- 
sumption when searching for a signal very close to the noise level. It searches for multiple 
transits spanning a large range of periods and transit start times. It was found from Monte 
Carlo simulations of our transit detection algorithm that the data must be analysed in two 
regions of the CMD, as they both have slightly different detection criteria. These criteria 
are described fully in Weldrake & Sackett (2004b), which deals with describing the MFA in 
detail. The resultant two bins are defined on Fig.l. Depths and durations of model transits 
are these expected for Hot Jupiters accounting for different positions the transit could have 
on the stellar disk as shown in Fig.l. Detailed testing revealed slightly different model depths 
within this range did not change the transit recoverability level to any significant degree, as 
the degree of scatter in the data is comparable to the transit depth. 

A total of 1,133,255 different transit models were compared to each sampled hghtcurve 
and a cross-correlation function (CCF) value created for each. A detection is indicated by 
a strong spike in the histogram of all CCF values. The results of the transit recoverability 
tests show that the code is very efficient at detecting transit signatures at all sampled orbital 
periods. For stars with photometric uncertainties <0.02 mag (rms), the recoverability spans 
95% to 35% for orbital periods of l-16d; including the recoverability of stars with rms 
uncertainty of 0.02-0.04 mag yields a weighted mean recoverability of 85% to 25% for the 
complete sample of stars over the same range of period. 



3.2. Expected Number of Hot Jupiter Detections 

Using the recoverability function R of the Monte Carlo simulations, seen in Fig. 11 of 
Weldrake & Sackett (2004b), we estimate the number of Hot Jupiter planets that should be 
detectable in our dataset, allowing a statistical significance to be placed on any null result. In 
order to do this, a number of factors need to be taken into consideration, including intrinsic 
frequency of Hot Jupiters, the total number of stars analysed, the geometric transit proba- 
bility at any particular given orbital period and the recoverability at that period delivered 
by the Matched Filter Algorithm. 

We begin by assuming that the frequency of Hot Jupiters, Ohj, in 47 Tuc is the same 
as in the solar neighbourhood, which has been estimated to be 0.8% (Butler et al. 2000; 
Basri 2004). This is then multiplied by the total number of main sequence stars available 
for analysis, n*, (21,920, 8,080 with rms <0.02 mag, and 13,840 with 0.02<rms<0.04 mag) 
to give the total number of Hot Jupiters that should be present in our dataset. The total 
number of Hot Jupiters are not randomly distributed in orbital period. Using the period 
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distribution^ of Hot Jupiters to date we split them into period bins of 0.5d width to separately 
calculate the expected number of detections in each period bin i and then sum these to arrive 
at: 



NhJ = Yl ^Hj{Pi) = (n. X Ohj) X [{fHj)i X {Ptran)i X {Ri)] (5) 
1=1 

where Nhj is the total number of detectable Hot Jupiters (HJ) present in the dataset, 
(n* X Ohj) is the total number of Hot Jupiters assumed to be present in the dataset, fnii 
is an estimate of the fraction of total Hot Jupiters expected to have periods in bin i as 
estimated from the solar neighbourhood, (Ptrani) is the geometric transit probability for bin 
i (Eq.4) and (Ri) is the recoverability delivered by our detection algorithm as determined in 
the Monte Carlo simulations, as seen in Fig. 11 of Weldrake & Sackett (2004b). 

The resulting expected distributions are presented in Table 2 and Fig.3. In these calcu- 
lations, the recent OGLE transit detections (Konacki et al 2003; Bouchy et al. 2004; Konacki 
et al. 2004) have been omitted, as the follow-up classification studies are on going, so the 
statistics cannot be calculated. Moreover, the OGLE detections were in environments other 
than the Solar Neighbourhood, and thus may have a different intrinsic Hot Jupiter frequency. 
Fig.3 shows that there should be a total of 7 detectable Hot Jupiters in our dataset, assuming 
that the formation probabihty is the same in 47 Tuc as in the solar neighbourhood. Using 
Poission statistics, a probability of no detections when 7 are expected would be 1/1,097, 
equating to a 3.3a null significance. 

It is noted that from recent transit discoveries this assumed radius may be anomalously 
large, and therefore was altered to a radius of 1.0 Rj. The resultant transit recoverability 
and expected numbers were not changed to any significant degree. 

Planet formation or survivability appears to depend on the metallicity of the host star 
(Gonzalez 1997; Laughlin 2000; Santos, Israehan, & Mayor 2001; Fischer & Valenti 2003), 
with the vast majority of planet host stars being of supersolar metallicity. However, only a 
few percent of the stars sampled for Hot Jupiter planets have metallicities similar to 47 Tuc 
(Fischer & Valenti 2003). Assuming that a better proxy for the frequency of Hot Jupiters 
in 47 Tuc is given by frequency of such planets around solar neighbourhood stars of similar 
metallicity, we show with the dotted line in Fig.3 that the expected number of planets in the 
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dataset would have to be revised downward to 0.9. Hence, a null result would be consistent 
with the expected number if the system metallicity is the dominant factor inhibiting planet 
formation in this environment. 



4. Results 

We ran our transit detection algorithm on the entire dataset of main sequence stars, 
with photometric rms <0.04 mag using slightly different detection criteria for the two main 
analysis bins outlined in Fig.l. For the 8,080 stars with rms <0.02 mag, a total of 64 possible 
"detections" were found. These detections were then phase-wrapped to all the significant 
periods to determine the nature of the periodicity that was seen. Phase wrapping indicated 
that the vast majority showed no clear pattern, but were instead caused by localised increases 
in scatter. These small-scale increases in random errors are very difficult to account for in 
the code when randomly spread in the timeseries, but are easily identified by eye when the 
data are phase wrapped to the periods singled out by the cross correlation function. 

In order for a candidate to be counted as a real transit detection, the lightcurve must 
pass several stringent criteria. The star must be within the V and V-I limits of the search. 
It must not be located close to a bright saturated object, as this can cause spurious transit 
signals in the data. It must also not be located very close to a bad column or other CCD 
blemishes, as image offsets during observing could cause the star to be periodically placed 
over these features. The star must have no heavily blended companions that degrade the 
ability of the photometric code to fit a reasonable point spread function (PSF) to the star. 
This can be tested easily by checking the goodness of fit for the PSF. A significantly high 
value indicates that the PSF is not being modelled appropriately, resulting in systematic 
deficiencies in lightcurve photometry. This could be caused by differing colours between the 
two heavily blended stars, or by other systematic effects introduced by ground-based pho- 
tometry that are difficult for the photometry code to handle. This was only a problem on 
very few of those stars with heavily blended neighbours. The candidates were comprehen- 
sively searched and those which displayed apparent transit features checked against these 
criteria. They were also checked with multiples of the detected periods, to see if any of the 
transits repeated. No signals which could be interpreted as being caused by a planet transit 
were observed. 

Our Matched Filter Algorithm did, however, detect some interesting systematic effects, 
including a star that was contaminated by a nearby detached eclipsing binary and gave the 
highest significance detection (9.5a) in this sample. The responsible binary was already 
identified as V39 (P=4.6015d) and published in WSBFa. The contaminated star is shown in 
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Fig.4, and is a good example of the ability of the code to detect low amphtude signals. The 
main eclipse occurs at Phase ($)=0 and $=1, with the depth and duration marked. The 
secondary eclipse can be clearly seen at $=0.5. No additional variable stars were identified in 
our transit search over the variables already identified in WSBFa using another algorithm. 
As only two binary stars were recovered, it seems that the vast majority of the variables 
available in this dataset (100 in total, 69 new discoveries) have already been identified and 
published in our companion paper WSBFa. 

Twelve transit-like features were identified in our search that had a common period of 
1.66d; one of these is plotted in Fig.5. This feature (at $=1) has the correct depth for a 
Hot Jupiter companion, and a duration (1.39 hours) that would be typical of a planet which 
only grazes the stellar disk. If this detection had not occured at a common periodicity, and 
at common places in the time-series, this would be a candidate for a grazing Hot Jupiter or 
eclipsing binary. We present this false detection here as an excellent example of the ability 
of the detection algorithm to identify even short duration transit-like events, which is typical 
of worst-case scenario transits where the planetary disk just fits inside the stellar disk. 

For the second sample of stars, with larger rms photometric uncertainties of 0.02-0.04 
mag, a total of 48 "detections" were found in the dataset. Again, these candidates were 

phase-wrapped at their respective determined periods and the nature of the periodicity seen. 
After checking in detail for the systematic effects listed above, no transits were observed. 

4.1. Interpretation 

Our null result has high statistical significance {3.3a), and is in agreement with the 
findings of Gilliland et al. (2000) that the environments of globular clusters (at least 47 
Tuc) lead to a dramatic reduction in Hot Jupiter formation probabilities. There are several 
possible scenarios that can explain this result. 

The main conclusion from the work of GiUiland et al. (2000) in the core of 47 Tuc was 
that either metallicity or crowding is the most important factor for significantly reducing the 
formation of Hot Jupiter planets in this environment. As our search is most sensitive to the 
uncrowdcd outer regions of the cluster, where the cluster relaxation time is longer, we would 
expect crowding to have a much smaller effect on our sample than samples taken from the 
cluster core. In any event, the work of Sigurdsson (1992) shows that the survivability of a 
planet at lAU semi-major axis is suprisingly long (~10^yrs) even for conditions typical of 
the core of 47 Tuc, and for P<5d orbital period planets, the stellar crowding would disrupt 
and eject only a small fraction from their parent systems. 
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Bonnell, Smith, Davies, & Horne (2001) undertook a further investigation into the 
effect of stellar crowding on planet formation. The conclusion was that the disruption of 
Hot Jupiter planets is unlikely to occur owing to encounters in a stellar cluster, due to 
the long timescales involved. Planets with long orbital period can however be ejected from 
their parent system. In the core of such a cluster with densities >10"^ star pc^^, even those 
planetary systems with semi-major axes of 0.1 au are likely to be disrupted and only those 
systems with semi-major axes <0.01 au are likely to survive. Lower densities typical of the 
halo of the cluster (where this work is most sensitive) will leave wider systems (<10 au 
semi-major axis) intact. 

Additionally, two of the variable stars (V78 and V93) identified in the field of 47 Tuc 
appear to be orbited by M-dwarfs with orbital periods of a few days (WSBFa). If these 
companions are indeed M-dwarfs this would further add to the argument that low mass 
objects can survive in the outer regions of 47 Tuc for appreciable times. 

Our null result therefore suggests that system metallicity is the dominant factor in- 
hibiting planet formation in 47 Tuc, and is consistent with the expected number of detected 
planets around stars of cluster metaUicity from radial velocity searches (Fig. 3). However, 47 
Tuc is an old system, and Vidal-Madjar et al. (2003) have shown that the transiting planet 
HD209458b seems to be evaporating its upper atmosphere to a significant degree. Could it be 
that all the planets have evaporated? These authors conclude that the evaporation timcscale 
is comparable to the main sequence lifetime of the host star, and these planets would then 
become detectable as hydrogen depleted objects. Lecavelier des Etangs ct al. (2004) give a 
range of evaporation rates for HD209458b. Even the highest of these (5x lO^^g/s^^) gives an 
evaporation timescale of ~8xl0^°yrs (assuming the rate remains constant), which is ~5.9 
times the age of 47 Tuc. 

To display our null result graphically, Fig. 6 shows the statistical significance of the result 
for a variety of different Hot Jupiter occurrence rates (fraction of stars with Hot Jupiter 
planets), plotted against orbital period. The currently accepted Hot Jupiter occurrence rate 
for the Solar Neighbourhood (Butler et al. 2000; Basri 2004) is plotted as a dotted line. We 
can hence rule out Hot Jupiter planets in 47 Tuc to 4(7 significance to P=3.6d, 3cr significance 
to a period of 5.5 days and 2a to P=9.2d, if the cluster has the same occurrence rate as the 
solar neighbourhood. From Fig.6, any occurrence rate hypothesis from 0.001 to 0.1 can be 
ruled out in 47 Tuc to a significance of la—>5a to a variety of orbital period limits. 
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5. Summary and Conclusions 

In this paper, we have presented results of a wide-field ground based search for Hot 
Jupiter transits in the globular cluster 47 Tucanae. The dataset in total included time 
series lightcurves (produced via difference imaging) of 21,920 main sequence stars of suitable 
photometric accuracy to search for planetary companions. The sampled orbital period range 
was l-16d, typical of Hot Jupiters found in the solar neighbourhood and defined by the 
length of our observing window. Seven planets were expected in the dataset, if the Hot 
Jupiter formation probability is the same as that in the solar neighbourhood, but none were 
found. This is a null result of high statistical significance (3.3(t), as tested with extensive 
Monte Carlo simulations on our dataset and is in agreement with the results of Gilliland et 
al. (2000) in an independant and complementary HST study of a region centered on the core 
of the cluster. 

The main aim of this work was to discriminate between stellar crowding and system 
metallicity as driving factors for the formation (or otherwise) of Hot Jupiter planets. By 
observing the whole visible extent of the cluster, (approaching 60% of the cluster tidal 
radius) and hence being more sensitive to regions of significantly less crowding than that 
of Gilliland et al. (2000), our results suggest that system metaUicity is the dominant factor 
determining the formation of Hot Jupiter planets in this environment. This result confirms 
the work of Gonzalez (1997); Laughhn (2000) and Santos, Israehan, & Mayor (2001), who 
note that extrasolar planets are preferrentially associated with stars of supersolar metallicity, 
and shows that planetary formation (at least for Hot Jupiters) is not a universal process, 
being strongly dependent on environment. 
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Fig. 1. — The expected depth and duration of planetary transits in 47 Tucanae involving 
main sequence stars. The predicted depths assume the planet passes on the 'average' chord 
length of the stellar disk, and has a period of 3.3d, with a radius of 1.3R/. The target 
range for this project is V=17.0-19.5, and hence has corresponding transit depths of 0.007 to 
0.031 magnitudes. The two defined boxes show the position of stars with rms photometric 
uncertainty <0.02 mag (top), and stars rms <0.04 mag (bottom), and were the two regions 
searched with the Matched Filter Algorithm. The uncertainly in our CMD dataset is also 



plotted as error bars. 
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Fig. 2. — The geometrical transit probability, Ptran for a solar radius primary in 47 Tuc, 
located in the middle of our transit search range, orbited by a 1.3Rj radius secondary. The 
probability decreases as period increases, and has been plotted to an upper limit of 16d, as a 
planet of this period could display three detectable (and hence confirmed) transits inside our 
observing window. By altering the assumed radius to IRj the resultant transit recoverability 
levels and expectations were not changed by any significant degree. 
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Fig. 3. — The number of expected Hot Jupiter planets detectable in our dataset, as a 
function of orbital period. These numbers incorporate the Monte Carlo simulation transit 
recoverability results, and are scaled with the distribution of Hot Jupiter planets in the Solar 
Neighbourhood. The dashed line is the expected number of planets for stars in the solar 
neighbourhood at the metaUicity of 47 Tuc, totaUing 0.9 planets. 
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Fig. 4. — The highest significance detection in the transit search for stars <0.02 mag rms. 
This is contamination of a neighbouring detached echpsing binary imposed on a 47 Tuc main 
sequence star, and is a good example of the abihty of the MFA to detect transit-hke features. 
The main eclipse occurs at $=0 and 1. The secondary eclipse can be seen at $=0.5. The 
photometric errorbars have been overplotted. 
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Fig. 5. — A typical detection which occurs at a common periodicity and common position 
in the time-series (at $=1) in the <0.02 mag rms sample. This detection has the expected 
depth of a Hot Jupiter companion, but only ~0.5 the duration. This example shows that 
a similarly deep transit but with twice the duration would be easily visible in the data. As 
this derived orbital period was common in the dataset, along with the MJD position of the 
'transit' features, this has been classified as a systematic effect. 
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Fig. 6. — To highlight the null result in a general way, this figure shows the statistical 
significance that can be placed on a range of Hot Jupiter occurrence rates (fraction of stars 
with Hot Jupiters) plotted against orbital period for 47 Tuc. The Hot Jupiter occurrence 
rate in the solar neighbourhood is plotted as a dotted line. Any Hot Jupiter occurrence rate 
hypothesis can be ruled out to the significance marked to the respective orbital period limit. 
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The parameters of 47 Tuc stars suitable for the transit search, as derived from the 
isochrones and luminosity functions of Bergbusch & VandenBerg (2001), kindly provided 
by R.Gilliland (private communication). The assumed 47 Tuc distance modulus is 13.4, as 

derived by those authors. 
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The expected number of detectable Hot Jupiter planets in the dataset as a function of 
orbital period P, Nhj{P), as described in section 3.2. The total number of expected 
detctions is 7 planets, assuming the formation probability of Hot Jupiters is the same in 47 

Tuc as it is in the solar neighbourhood. 



